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Abstract 

It is argued that the observed dehcit of solar and atmospheric neutrinos can be ex¬ 
plained by neutrino oscillations ^ and —>■ z/' involving two hypothetic sterile 

neutrinos z/^ and z/' (blind to all Standard-Model interactions). They are keen to mix 
nearly maximally with z/g and z/^, respectively, to form neutrino mass states z/i, z /4 and 
Our argument is presented in the framework of a model of fermion ’’texture” for¬ 
mulated previously, which implies the existence of two sterile neutrinos beside the three 
conventional. 
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1. Introduction 


The recent findings [1] of Super-Kamiokande atmospheric-neutrino experiment 
brought to us the important message that the observed deficit of atmospheric z/^’s seems 
to be really caused by neutrino oscillations, related to nearly maximal mixing of with 
another neutrino. This may be Vr or, alternatively, a new sterile neutrino (blind to 
all Standard-Model interactions). The z/g neutrino is here excluded from being a mixing 
partner of z/^ by the negative result of CHOOZ long-baseline reactor experiment [2] which 
found no evidence for the disappearance modes of Vf., in particular Vf. ^ parameter 

region overlapping the range of sin^ 26*atm and Amatm observed in the Super-Kamiokande 
experiment. 

The survival probability for z/^, when analized experimentally in two-flavor form 

P ^ z/^) = 1 - sin^ 26»atm sin^ (l.27Amat„i L/E^ , (1) 

leads to the parameters [1] 

sin2 20atm = 0(1) ~ 0.82 to 1 (2) 

and 

Am^^m ~ (0-5 to 6) X 10“=^ eV^ (3) 

at the 90% confidence level (note that the value Amatm ~ 5 x 10“^ eV^ corresponds to 
the lower limit of the previous Kamiokande estimate of [3]). If z/^ is responsible 

for this nearly maximal mixing of z/^, then the disappearance probability for z/^ in the 
mode z/^ is 


p i^r) = sin^ 26'atm slu^ (l.27Amatm L/e) . (4) 

In the present paper, we conjecture that it is rather a sterile neutrino (denoted here 
by z/() which is responsible for such a nearly maximal mixing of z/^ (whet her it is not 
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or is Vr constitutes a crucial point of our conjecture which, unfortunately, is not at the 
moment easy to decide experimentally [1]). We conjecture moreover that another sterile 
neutrino (denoted by Vs) mixes nearly maximally with z/g, causing the observed deficit 
of solar z/e’s. In such a way, we introduce a unified picture of neutrino oscillations as 
being related to nearly maximal mixing of two sterile neutrinos Vs and z/( with z/g and 
z/^, respectively. Of course, this mixing of z/^ and z/( is not forbidden by the weak isospin 
Ja and weak hypercharge Y of Vf. and z/^ , as the conservation of these weak charges is 
spontaneously broken, except for their combination Q = + Y/2 (equal to zero for z/g 

and z/^). We should like also to remark that the sterile neutrinos Vg and z/(, interacting 
only gravitionally, would be responsible for the existence of a Standard-Model-inactive 
fraction of the dark matter. 

Note that the existence of just two sterile neutrinos (blind to all Standard-Model 
interactions), beside three families of Standard-Model-active leptons and quarks, turns 
out to be natural in the model of lepton and quark ’’texture” we develop since some time 
[4,5] (c/. Eqs. (A.15) in Appendix). In this model, all neutrinos are Dirac particles having 
both lefthanded and righthanded parts. 

For the Standard-Model-active neutrinos z/g, z/^, z/,-, charged leptons e“ , /i“ , r“, 
up quarks u , c, t and down quarks d, s , b we came to a proposal [5] (c/. Eq. (A.10) in 
Appendix) of unified algebraic structure of their mass matrices (^M-p^ (/ = z/ , e , u , d) 
in the three-dimensional family space (i, j = 1, 2 , 3). In the case of leptons (/ = z/, e), 
this proposal reads 




0 \ 



4/i(A(80 + e(/)2)/9 




(5) 


\ 0 24/i(A(624 + e(A2)/25 j 

Here, and denote real constants to be determined from the present 

and future experimental data for lepton masses and mixing parameters and are 
mass-dimensional). 

For charged leptons, when assuming that the off-diagonal elements of the mass matrix 
pMpP given in Eq. (5) can be treated as a small perturbation of its diagonal terms, we 
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calculate in the lowest (quadratic) perturbative order in [5]: 
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When the experimental nie and [6] are used as inputs, Eqs. (6) give [5] 


nir 


ye) 2 



1776.80 + 10.2112 



MeV , 


0.172329 + 0 



85.9924 MeV + O 



(7) 


We can see that the predicted value of rrir agrees very well with its experimental hgure 

^exp ^ 

MeV [6], even in the zero-order perturbative calculation. To estimate 

(^a ie)/ye)y 

, we take this experimental figure as another input. Then, 

(y^) “ “'“OiaSi . (8) 

SO it is not inconsistent with zero. 

The unitary matrix diagonalizing the mass matrix according to the 

relation M^) = diag(rrie, nir), assumes in the lowest (quadratic) pertur¬ 
bative order in the form 
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where the small is neglected. Of course, in the limit of 0, we obtain (U-f^ 

{S^Y■ 

For neutrinos, we will assume in this paper that is very small and 


= 0 , ( 10 ) 
in contrast to the possibility of Y 0 charged leptons [cf. Eq. ( 8 )]. Then, for 
conventional neutrinos = (5p) and so, Ug, z/^, can mix only by means of the 

trivial lepton CKM matrix (Vij) = what is a minor effect, 

vanishing in the limit of 0. Instead, allowing in this paper for the existence of 

two sterile neutrinos Ug and z/(, we will extend the 3x3 neutrino mass matrix 
{i,j = 1,2,3), given through Eqs. (5) and (10), to a 5 x 5 neutrino mass matrix 
(/, J = 1, 2, 3,4, 5) with = Mj'j*. Explicitly, we will assume that 
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( 11 ) 


and ~ /7, ~ consequence of Eqs. (A.19) and (A.20) {cf. 

Appendix). It will turn out that the matrix elements = M^* and = M^* 
lead to the mixing of neutrino flavor states Ve with Vg and with within neutrino mass 
states z/i, z /4 and z/ 2 , r's, respectively. 


2. Neutrino mass states 

The eigenvalues of the extended mass matrix given in Eq. (11) are Dirac 

masses of hve neutrino mass states z/i, z /2 , ^3 , ^4 , They are 
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In Section 4, the masses and rriy^ will tnrn ont to be negative, what is irrelevant in 
the case of Dirac particles for which only masses squared are measurable (so, and 

\mu^\ will be the phenomenological masses of and z/ 2 ). 

The corresponding 5x5 unitary matrix diagonalizing the mass matrix 

according to the equality = diag{m^-^ , , takes the 

form 




The neutrino flavor states = Ug, , Ug, (of which z/g, z/^ , Vti or rather 

their lefthanded parts, stand for the observed weak-interaction neutrino states and Vg, 
v'g denote their unobserved sterile partners) are related to the neutrino mass states z// = 
5 ^^2 , r'a , ^4 , r's through a hve-dimensional unitary transformation 


j 


( 15 ) 
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with {y*j^ = Here, 
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where (Uif^ is the charged-lepton diagonalizing matrix given in Eq. (9) and 


t/h) _ 0 = 


^ 24 


25 


u^} = 0 = 


"4j 


5j 


/•/(e) _ 1 _ Trie) 

O 44 — 1 — O 55 


( 17 ) 


The last equations follow from the fact that charged leptons get no sterile partners. Thus, 
from Eq. (16) 


T/ _ TtM *TTi<^) T/ _ TjM* T/ _TtW)* 
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( 18 ) 


Of course, the 5x5 unitary matrix {Vaj) is a five-dimensional lepton counterpart of the 
familiar CKM matrix for quarks. The charged leptons e~ , ii~ , t~ are here counterparts 
of the up quarks u, c,t (both with diagonalized mass matrix). 

From Eqs. (18), with the use of Eqs. (13) and (9), we can calculate the matrix elements 
Vaj in the lowest (quadratic) perturbative order in Q;h)/^(®). Writing for convenience 
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and 
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Vu = 

V 24 = 0 
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V"i5 = 0 
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In the limit of 0, the only nonzero matrix elements V^j are 
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3. Neutrino oscillations 

Having once found the elements (19) and (20) of the extended lepton CKM matrix, we 
are able to calculate the probabilities of neutrino oscillations (in the vacuum), 

using the familiar formula: 

where //^(O) = z/q, {up\ = (0|z/^ and {i'p\va) = ^/ 3 q- Here, as usual, f/|p| = LjE (c = 
1 = h), what is equal to 4 x 1.2663L/i? if — ml^ , L and E are measured in eV^, 
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m and MeV, respectively. Of course, L is the source-detector distance (the baseline). In 
the following, it will be convenient to denote 


xlk = 1.2663 
and use the identity cos 2 xlk = 1 — 2 sin^ xlk 




E 


(24) 


From Eqs. (23), (19) and (20) we derive by explicit calculations the following neutrino- 
oscillation formulae valid in the lowest (quadratic) perturbative order in 
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In the limit of 0, the only nonzero neutrino-oscillation probabilities are 


(25) 


y2 

P {l^e ^ l^s) ^ 4 ^^ ^^ 2)2 ^41 , 

4 (1 ^^ 2)2 sin^^52 • (26) 

The formulae (25) for the disappearance modes of Ug and u^ imply the following 
survival probabilities for u^ and z/^: 



P(z/e ^ l^e) = ^-P{^e ^ ^ I^t) - P (l^e ^ Vs)-P{Ve ^ Pg) 
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In the limit of 0, we obtain 


P (z/e ^ Z/e) ^ 1 - 4 


y2 


(1 + F 2)2 
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(29) 


and 


X2 

p ^ 12/.) ^ 1 - 4 ^^ ^^2)2 ^52 • (30) 

The last two formulae are to be compared with solar-neutrino and atmospheric- 
neutrino experiments, respectively. 

4. Atmospheric and solar neutrinos 

In the case of atmospheric neutrinos, we compare our formula (30) with Eq. (1). 
Then, for instance. 


(more generally: ~ 0.82 to 1) and 


4X2 

(1+X2)2 


~ 0.9 


(31) 
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(32) 


ml^ — ml^ ~ 5 X 10 ^ eV^ 

(more generally (0.5 to 6) X 10-3 eV^). 

From the input (31) we get 


X ~ 0.721 

and, through the second Eq. (14), 

_ 1 - X2 1 
2|m£)| “ 2X ^3 

or 
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On the other hand, the third mass formula (12) and the input (32) give 


(33) 


(34) 


(35) 


(a4"> + a 4') \/(A4>-A4>)" + 4|A4y = ml - m^ ~ 5 x IQ-* eV (36) 
or, with the use of Eqs. (34) and (33), 


{rnl^ - rnl^) ~ 1-58 x 10 ^ eV^ . 

With the formulae ~ 320/i*^^V261 and ~ 48/i*''^V'^ 'w® have 
45.5/i‘^'^^^. Hence, Eq. (37) leads to 
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Then, 


Mg) ~ 7.25 X 10-3 eV , Mg) ~ 4.04 x 10-^ eV (39) 

and so, from Eq. (35) 


I Mg) I ~ 4.97 X 10-2 eV. 


(40) 
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Finally, with the values (39) and (40) the third mass formula (12) gives 

-2.86 X 10-2 eV 
7.62 X 10-2 eV ‘ 

In this way, all parameters appearing in our model of neutrino ’’texture”, needed to 
explain the observed dehcit of atmospheric z/^’s in terms of neutrino oscillations z/(, 
are determined. 

In the case of solar neutrinos, we compare our formula (29) with the survival proba¬ 
bility for z/e, usually analized experimentally in two-flavor form 

P(z/e —z/e) = 1 — sin2 26*soisin2 ^1.27AmgoiL/i7^ . (42) 

Taking into account the so-called vacuum fit [7] (he., one that is not enhanced by the 
resonant MSW mechanism [8] in the Sun matter), we have the parameters 




sin2 26*soi ~ 0.65 to 1 , ~ (5 to 8) x 10 eV^ , (43) 

what shows a large mixing and a very small difference of masses squared. Then, for 
instance. 


(more generally: ~ 0.65 to 1) and 
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From the input (44) we obtain 
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On the other hand, the first mass formula (12) and the input (45) lead to 


(A/if + mIJ') -a4;>)04|mS'|2 = ml^ - mj, ~ 7 x IQ-" eV'^ 


(49) 


or, through Eqs. (47) and (46), to 
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Eventually, with the values (52) and (53) the hrst mass formula (12) implies 




(53) 


(54) 


-1.32 X 10-® eV 
8.48 X 10“® eV ' 

V 

In such a way, all parameters contained in our model of neutrino ’’texture”, needed to 
describe the observed dehcit of solar u^s in terms of neutrino oscillations Vg in the 
vacuum, are determined. 

Our last item is concerned with the LSND accelerator experiment that reported the 
detection of —> Ve and Ve oscillations by observing z/g’s and z/g’s in a beam 

of P^’s and z/^’s produced in ■n~ and tt’*' decays, respectively [9]. The observed excess 
of Pg’s and z/g’s, analized in terms of two-flavor neutrino-oscillation formula, implies a 
considerable amplitude sin^26*LSND, too large to be explained by our formula (25) for 


P(z/^ —^ ^e) = P{^e where the amplitude at sin^a; 2 i. 


16 
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is small; 
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as it follows from Eq. (8). Here, the central value is 
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I would like to thank Jan Krolikowski for several helpful discussions. 


Appendix: Unified ’’textnre dynamics” 


In this Appendix the idea of a model of fermion ’’texture” that we develop since some 
time [4,5] is outlined. In particular, the existence of two sterile neutrinos Vg and z/' turns 
out to follow naturally. 

Let us introduce the following 3x3 matrices in the space of three fermion families: 
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With the matrix 


n = a^a 



they satisfy the commutation relations 


a, h] = a , [a^, n] = —a^ 


(A.l) 


(A.2) 


(A.3) 


characteristic for annihilation and creation matrices, while n plays the role of an occup¬ 
ation-number matrix. However, in addition, they obey the ’’truncation” identities 


a3 = o,a^3 = 0. (A.4) 

Note that due to Eqs. (A.4) the bosonic canonical commutation relation [a, al] = 1 does 
not hold, being replaced by the relation [a, al] = diag (1,1, —2). 
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In consequence of Eqs. (A.l), (A.2) and (A.3), we get n\n) = n\n) as well as a\n) = 
^Jn\n — 1) and a)\n) = \/n + l\n +1) (n = 0,1, 2), however, al|2) = 0 (ie., |3) = 0) in 
addition to a^|0) = 0 (ie., | — 1) = 0). Evidently, n = 0,1, 2 may play the role of a vector 
index in our three-dimensional matrix calculus. 

It is natural to expect that the Gell-Mann matrices (generating the horizontal SU(3) 
algebra) can be built up from a and ah In fact. 
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(A.5) 


Inversely, a = (Ai + iA 2)/2 + \/2(\e + *7)/2 and a) = (Ai — 1 X 2 )/'^ + v^(\ — iXj)/2. 
A message we get from these relationships is that a horizontal field formalism, always 
simple (linear) in terms of A^ (A = 1, 2,..., 8) and 1, is generally not simple in terms of 
a and ah In particular, a nontrivial SU(3)-symmetric horizontal formalism is not simple 
in a and ah Inversely, a nontrivial horizontal held formalism, if simple (linear and/or 
quadratic and/or cubic) in terms of a and a\ cannot be SU(3)-symmetric. 

Now, let us consider the following ansatz [5]: 


MU) = p^A/)(/)^i/2 (^f = u , e , u , d) , 


(A.6) 
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where 
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It is the matter of an easy calculation to show that the matrices (A.6) get explicitly the 
form [5]: 



1 p(7)£(/)2 


0 ^ 

M)7) = 7 

29 


477)(80 + e)7)2)/9 

873(a)7)_/3(/))e*^‘^^ 


1 0 
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In this paper we write also = (^M-p^ {i,j = 1, 2, 3). 

In a more detailed construction following from our idea about the origin of three 
fermion families [4], each eigenvalue A^ = 1, 3 , 5 of the matrix N corresponds (for any / = 
u , e , u , d) to a wave function carrying A^ = 1, 3 , 5 Dirac bispinor indices: ai, 0 : 2 ,..., av 
of which one, say ai, is coupled to the external Standard-Model gauge helds, while the 
remaining A^ — 1 = 0, 2,4: q; 2 ,---,« 7 v (that are not coupled to these helds) are fully 
antisymmetric under permutations. So, the latter obey Fermi statistics along with the 
Pauli principle implying that really N — 1 <4, because each ctj = 1,2, 3,4. Then, the 


16 



three wave functions corresponding to TV = 1,3,5 can be reduced to three other wave 
functions carrying only one Dirac bispinor index ai (and so, spin 1/2), 


i’ll, = i’al , 

'f) = —e 

oi T ttlQ:2a3a4a5 

and appearing (up to the sign) with the multiplicities 1, 4 and 24, respectively. In this 
argument, for N = 3 the requirement of relativistic covariance of the wave function (and 
the related probability current) is applied explicitly [4]. The weighting matrix as 
given in Eq. (A. 7) gets as its elements the square roots of these multiplicities, normalized 
in such a way that Tr p = 1. 

In Eqs. (A.11), the indices a, (i = 1,2,... ,A^) are of Jacobi type: ai is a ’’centre- 
of-mass” Dirac bispinor index, while 0 : 2 ,..., av are ’’relative” Dirac bispinor indices. 
In fact, ai {i = 1,2,..., N) are dehned by chiral representations of Tf matrices (i = 
1,2,..., N) being the (properly normalized) Jacobi combinations of some individual 7 / 
matrices {i = 1,2,..., N), where, in particular, T^ = (1/V^) [4]. For them 

|r/,r/| = 2Sijg’"' {i,j = 1,2,..., N), in consequence of the anticommutation rela¬ 
tions { 7 *^, 7 /} = valid for any 7 / and 7 J. Then, the Dirac-type equations 

(Ti ■ [p — gA{x)] — M} ijj{x) = 0 {N = 1,2,3,...) [4], independent of T^ ,..., Tp^, hold for 
the fundamental-particle wave functions ipi^x) = (t/’aio 2 ...ajv(^))) where A^ = 1,3,5 in the 
case of fermion wave functions (A.11). Here, pTi ■ A{x) symbolizes the Standard-Model 
coupling. 

Note that all four matrices (^f = u, e, u, d) dehned by Eqs. (A.6) — (A.9) 

and (A.l) have a common structure, differing from each other only by the values of their 
parameters aA\ (3A\ CA'’ and We proposed the fermion mass matrices 

to be of this unihed form [5]. Then, Eqs. (A.6) and (A.8) dehne a quantum-mechanical 
model for the ’’textnre” of fermion mass matrices (/ = u, e,u, d). Such an 

approach may be called ’’textnre dynamics”. 

The fermion mass matrix containing the kernel h^A given in Eq. (A.8), consists 


_,/,(/) _,/,(/) 

— Ya, 12 — Ya, 34 ) 

= i’aluSi ) (A. 11) 
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of a diagonal part proportional to and of an off-diagonal part involving linearly 
and The off-diagonal part of describes the mixing of three eigenvalues 


/i 


(/) 




+ (£<•'12 - l) N-'^ + (N = 1,3,5) 


(A.12) 


of its diagonal part. Beside the term that appears only for N = 5, each of these 

eigenvalues is the sum of two terms containing They are: (i) a term that 

may be interpreted as an ” interaction” of N elements (” intrinsic partons”) treated on the 
same footing, and (u) another term 


^(/) |^£(/)2 _ p2 ^ [N\/{N - 1)!]“^ = N-^ (A.13) 

that may describe an additional ’’interaction” with itself of one element arbitrarily chosen 
among N elements of which the remaining iV — 1 are undistinguishable. Therefore, 
the total ’’interaction” with itself of this (arbitrarily) distinguished ’’parton” is 
(£(/)2 _ i)_/V“^], so it becomes in the hrst fermion family. 

The form (A. 11) of three fermion wave functions shows that that each ’’intrinsic par- 
ton” carries a Dirac bispinor index (of the Jacobi type). For the (arbitrarily) distinguished 
’’parton”, this index, considered in the framework of a fermion wave equation, is coupled 
to the external gauge helds of the Standard Model. Thus, this ’’parton” carries the total 
spin 1/2 of the fermion as well as a set of its Standard-Model charges corresponding to 
f = u , e , u , d. For the A^ —1 undistinguishable ’’partons”, obeying Fermi statistics along 
with the Pauli principle, their Dirac bispinor indices are mutually coupled, resulting into 
Lorentz scalars, while their number A^ — 1 = 0, 2, 4 differentiates between three fermion 
families (for each f = u, e, u, d). These ’’partons” are free of Standard-Model charges. 

Evidently, the intriguing question arises, how to interpret two possible boson families 
corresponding to the number A^ — 1 = 1,3 of undistinguishable ’’partons” [10]. In the 
present paper this problem is not discussed. Here, we would like only to point out that 
three fermion families A^ = 1, 3, 5 differ from these two hypothetic boson families A^ = 2,4 
by the full pairing of their A^ — 1 = 0,2,4 undistinguishable ’’partons”. So, the boson 
families, containing an odd number N — 1 = 1, 3 of such ’’partons”, might be considerably 
heavier. Note that the wave functions corresponding to A^ = 2,4 can be reduced (under 
some relativistic requirements) to two other wave functions carrying only spin 0, 
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1 

1 


{c-V) 

'^0:10:20:30:4 


0102 




(/) 

0102 




if) 

O 1 O 2 O 3 O 4 


= ^ - 



(A. 14) 


and appearing (up to the sign) with the multiplicities 2 and 6, respectively. 

Another important question also appears, namely, what is the interpretation of two 
fermions corresponding to the number iV = 1, 3 of undistinguishable ’’partons” only. Such 
fermions can carry exclusively spin 1/2 (for N = 3: under some relativistic requirements). 
Of course, they are free of Standard-Model charges and so, can be considered as two sterile 
neutrinos with the wave functions 


^sai — 01 — f/oi ) 


— t/soi = X i^ '1 ) ^02030405^' 

^ h ^ / Qi 02 


02030405 "r 030405 


^/>134 for tti = 1 
-^234 for tti = 2 
■^312 for tti = 3 (A. 15) 

-'«/’4i2 for tti = 4 


appearing (up to the sign) with the multiplicities 1 and 6, respectively. 

For these sterile neutrinos one may introduce the 2x2 mass matrix = 

^3 )i/2))(s)^s)i/2 ^ where 

while the diagonal part of is conjectured to have the eigenvalues 




is) 


N‘ 




p2 


with Pn = N\/N\ = 1 (iV = 1, 3) 


(A.17) 


Now, one ’’intrinsic parton” is arbitrarily chosen (to carry the total spin 1/2 of the fermion) 
among N ’’intrinsic partons” that all are undistinguishable [in contrast to Eqs. (A. 12) 
and (A. 13)]. This gives the diagonal part of equal to 


1 f Q 

7!^ 0 + 


(A.18) 
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Thus, the diagonal matrix elements and of the 5x5 neutrino mass matrix 
(/, J = 1, 2, 3,4, 5) introduced in Eq. (11) get the forms 


Mi:> = 0 , a4' = ^ (8 + 


7 


■ 00 

with expected to be very small. In the present paper we will assume that 


(A.19) 




^is)2 ^^iu) 2 _ 


(A.20) 


in Eqs. (A.19). 

The possibility of existence of two bosons corresponding to the number N = 2,4 of 
undistinguishable ’’partons” only ought to be also considered. Such bosons can carry 
exclusively spin 0 (for N = 2: under some relativistic requirements). Obviously, they are 
free of Standard-Model charges and so, may be considered as two ’’sterile scalars” with 
the wave functions 


02 — . (C y ') '00102 '012 034 ) 04 — „ . ^«lQ2Q3Q40ciia2a3Ci4 01234 

4 \ / 0102 24 

appearing (up to the sign) with the multiplicities 4 and 24, respectively. 

A priori, the ’’intrinsic partons” may be either strictly algebraic objects providing 
fundamental fermions (leptons and quarks) with new family degrees of freedom, or may 
give us a signal of a new spatial substructure of fundamental fermions (built up of spatial 
’’intrinsic partons” = preons, related to the individual yf as well as xf and pf {i = 
1,2,...,N)] note that here yf ’s anticommute for different i !). Our idea about the 
origin of three fermion families [4] chooses the hrst option. The difficult problem of new 
non-Standard-Model forces, responsible for the binding of N preons within fundamental 
fermions, does not arise in this option. 

However, if the second option is true, then this irksome (though certainly profound) 
problem does arise and must be solved. 


20 



References 


1. Y. Fukuda et al. (Super-Kamiokande Collaboration) , ’’Evidence for oscillation of 
atmospheric neutrinos”, to appear in Phys. Rev. Lett.] and references therein. 

2. M. Appolonio et al. (CHOOZ Collaboration), Phys. Lett. B 420, 397 (1998). 

3. Y. Fukuda et al. (Kamiokande Collaboration), Phys. Lett. B 335, 237 (1994). 

4. W. Krolikowski, Acta Phys. Pol. B 21, 871 (1990); Phys. Rev. D 45, 3222 (1992); in 
Spinors, Twistors, Clifford Algebras and Quantum Deformations (Proc. 2nd Max Born 
Symposium 1992), eds. Z. Oziewicz et al, 1993, Kluwer Acad. Press. 

5. W. Krolikowski, Acta Phys. Pol B 27, 2121 (1996); B 28, 1643 (1997); B 29, 629 
(1998); B 29, 755 (1998); |hep-ph/9803'^ . 

6. Review of Particle Physics, Phys. Rev. D 54, 1 (1996), Part I. 

7. N. Hata and P. Langacker, |hep-ph/9705339| ; cf. also G.L. Fogli, E. Lisi and D. Mon- 
tanino, hep-ph /9709473. 

8. L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); S.P. Mikheyev and A. Y. Smirnow, Sov. 
J. Nucl. Phys. 42, 913 (1985); Nuovo Cimento, C 9, 17 (1986). 

9. C. Athanassopoulos et al (LSND Collaboration), Phys. Rev. C 54, 2685 (1996); Phys. 
Rev. Lett. 77, 3082 (1996); |nucl-ex/9709~00^ . 

10. W. Krolikowski, Phys. Rev. D 46, 5188 (1992); Acta Phys. Pol. B 24, 1149 (1993); 
B 26, 1217 (1995); Nuovo Cimento , 107 A, 69 (1994). 


21 









